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Rayleigh - Taylor Gravity Waves and Quasiperiodic Oscillation 

Phenomenon in X-ray Binaries 

Lev Titarchuk 1,2 

ABSTRACT 

Accretion onto compact objects in X-ray binaries (black hole, neutron star 
(NS), white dwarf) is characterized by non-uniform flow density profiles. Such 
an effect of heterogeneity in presence of gravitational forces and pressure gra- 
dients exhibits Raylegh-Taylor gravity waves (RTGW). They should be seen as 
quasioperiodic wave oscillations (QPO) of the accretion flow in the transition 
(boundary) layer between the Keplerian disk and the central object. In this pa- 
per I show that the main QPO frequency, which is very close to the Keplerian 
frequency, is split into separate frequencies (hybrid and low branch) under the 
influence of the gravitational forces in the rotational frame of reference. The 
RTGWs must be present and the related QPOs should be detected in any sys- 
tem where the gravity, buoynancy and Coriolis force effects cannot be excluded 
(even in the Earth and solar environments). The observed low and high QPO 
frequencies are an intrinsic signature of the RTGW. I elaborate the conditions 
for the density profile when the RTGW oscillations are stable. A comparison of 
the inferred QPO frequencies with QPO observations is presented. I find that 
hectohertz frequencies detected from NS binaries can be identified as the RTGW 
low branch frequencies. I also predict that an observer can see the double NS spin 
frequency during the NS long (super) burst events when the pressure gradients 
and buoyant forces are suppressed. The Coriolis force is the only force which 
acts in the rotational frame of reference and its presence causes perfect coherent 
pulsations with a frequency twice of the NS spin. The QPO observations of neu- 
tron binaries have established that the high QPO frequencies do not go beyond 
of the certain upper limit. I explain this observational effect as a result of the 
density profile inversions. Also I demonstrate that a particular problem of the 
gravity waves in the rotational frame of reference in the approximation of very 
small pressure gradients is reduced to the problem of the classical oscillator in 
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the rotational frame of reference which was previously introduced and applied 
for the interpretation of kHZ QPO observation by Osherovich & Titarchuk. 

Subject headings: Accretion, accretion disks — (magnetohydrodynamics:) MHD — 
stars:oscillations (including pulsations) — stars: neutron— X-ray:binaries 


1. Introduction 

The theory of oscillations of rotating fluids of variable density [Rayleigh- Taylor (R-T) 
effect] was developed in detail by Chandrasekhar (1961), hereafter C61. A large variety of 
the magnitohydrodynamic (MHD) problems, including the stability of inviscid and viscous 
fluids in the case of two uniform layers separated by a horizontal boundary with and without 
rotation as well as the magnetic field effects were analyzed using perturbation technique. 
The simplest case of a one-dimensional gravitational force was studied. A similar analysis 
was also implemented for the case of an exponentially varying density. It follows from C61 
that the quasiperiodic oscillations (QPO) with the twin “kilohertz” frequencies represent 
the main frequencies of the stable gravity waves in the rotational frame of reference. It 
also as follows from C61 that the twin “kilohertz” frequencies should be on the order of the 
Keplerian frequency. 

In this Paper I present the results of the study of the R-T effect for a particular case of 
fluid oscillations in an accretion flow under influence of a central gravitational force. This 
particular R-T analysis is important in view of the high and low frequency detection by the 
Rossi X-ray Timing Explorer (RXTE) in a number of low mass X-ray binaries (Strohmayer 
et al. 1996, van der Klis et al. 1996), black hole candidate sources (Morgan, Remillard & 
Greiner 1997; Strohmayer 2001 a, b; Remillard 2002) and by Extreme Ultraviolet Explorer, 
Chandra X-ray observatory and optical observations in white dwarfs (Mauche 2002 and 
Woudt k Warner 2002). The presence of two observed peaks with frequencies and u 2 
in the upper part of the power spectrum became a natural starting point in modeling the 
phenomena. In NS binaries, for example Sco X-l, the lower frequency part of the power 
spectrum, contains two horizontal branch oscillation (HBO) frequencies i'hbo ~ 45 Hz and 
V2HBO ~ 90 Hz (probably the second harmonic of i 'hbo) which slowly increase with the 
increase of rq and v 2 (van der Klis et al. 2000). Any plausible model faces the challenging 
task of describing the dependences of the peak separation A u = v 2 — v\ on iq and v 2 . 
Attempts have been made to relate u\ and u 2 and the peak separation A v = v 2 — v\ with 
the neutron star (NS) spin. In the sonic point beat frequency model by Miller, Lamb k 
Psaltis (1998) the kHz peak separation Av is considered to be close to the NS spin frequency 
and thus Av is predicted to be constant. However observations of kHz QPOs in a number 
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of binaries (Sco X-l, 4U 1728-34, 4U 1608-52, 4U 1702-429 and etc) show that the peak 
separation decreases systematically when the high [kilohertz (kHz)] frequencies increase (for 
a recent review see van der Klis 2000, hereafter VDK). For Sco X-l VDK found that the 
peak separation of kHz QPO frequencies changes from 320 Hz to 220 Hz when the lower kHz 
peak U\ changes from 500 Hz to 850 Hz. 

The correlation between high frequency (lower kHz frequency) and low frequency (broad 
noise component) QPOs previously found by Psaltis, Belloni &; van der Klis (1999) for 
black hole (BH) and neutron star (NS) systems has been recently extended over two orders 
of magnitude by Mauche (2002) to white dwarf (WD) binaries. Accepting the reasonable 
assumption that the same mechanism produces the QPO in WD, NS and BH binaries, one 
can argue that the data exclude relativistic models and, beat frequency models as well as any 
model requiring either the presence or absence of a stellar surface or a strong magnetic field. 

The transition layer model (TLM) was introduced by Titarchuk, Lapidus & Muslimov 
(1998), hereafter TLM98, to explain the dynamical adjustment of a Keplerian disk to the 
innermost sub-Keplerian boundary conditions (it is at the star surface for NSs and WDs). 
TLM98 argued that a shock should occur where the Keplerian disk adjusts to the sub- 
Keplerian flow. Thus the transition layer bounded between the sub-Keplerian boundary and 
the adjustment radius can undergo various type of oscillations under the influence of the 
gas, radiation, magnetic pressure and gravitational force. Osherovich & Titarchuk (1999), 
hereafter OT99, suggested that the phenomological model of a one-dimensional classical 
oscillator in the rotational frame of reference could explain the observed correlations between 
twin kHz frequencies iq, tq and the HBO frequencies. They further suggested that the 
oscillations of the fluid element that bounced from the disk shock region (at the adjustment 
radius) would be seen as two independent oscillations parallel and perpendicular to the disk 
plane respectively. This is due to the presence of a Coriolis force in the magnetospheric 
rotational frame of reference. In this paper I show that when the pressure gradients can 
be neglected, the problem of the Rayleigh-Taylor wave oscillations (buoynancy effect) in 
the rotational frame of reference is reduced to the OT99 formulation (see §3). This result 
provides a solid basis for application of the OT99 model for interpretation of the QPO 
phenomena observed in X-ray binaries. 

The main goal of this Paper is to demonstrate that there is an inevitable effect of 
the gravity wave oscillations in the heterogeneous fluid of the accretion flow near compact 
objects. In §2, I formulate the problem of the gravity wave propagation in the bounded 
medium in the rotational frame of reference. In §3 I present analysis and solutions of this 
problem using a perturbation method in the context of three dimensional periodic waves 
with assigned wave numbers. Applications of the gravity modes and their relations with the 


- 4 - 


QPO observations is presented in §4. Summary and conclusions are drawn in §5. 


2. The Problem of Gravity Wave Propagation in a Bounded Medium 


Lord Rayleigh (1883) has treated the non- rotating inviscid case of the present problem. 
He developed a general theory for any density configuration po(z) (the z axis being the 
upward drawn vertical). Rayleigh’s treatment of inviscid superposed fluids was extended 
by Bjerkness et al. (1933) to include the influence of rotation. With an assumption of the 
gravitational force directed vertically, Hide (1956) developed the theory for any density Po{z) 
and viscosity po{z) profiles and any <5 (where 8 is an angle between the rotational axis and 
the vertical). In this study we consider a case of an inviscid and incompressible fluid. 

The equation of relative motion appropriate to the problem is 


dui d 

p ~dt +pu 'd^ Ui 


n dp 

2pyle i j k UjS( C ~3x~- 


(i) 


where the fluid is supposed to rotate uniformly about an axis whose direction is specified by 
the unit vector s = Q/Q. The tensor notation follows the summation convention and the 
unit vector e in the radial direction is introduced. In equation (1) p denotes the density, fi 
the angular velocity of rotation, tq the ith component of the (Eulerian) velocity vector, p is 
the pressure and g is the acceleration due to gravity. This term g can represent the effective 
gravity which includes the centrifugal and radiation pressure forces. For an incompressible 
fluid the continuity equation is 


out 

d^l~ ' 


(2) 


Because (similar to C61)diffusion effects are ignored in this analysis, an individual fluid 
element retains the same density throughout its motion. Hence it is required that 


Dp _ dp dp _ 
Dt dt Uj dxj 


( 3 ) 


Because the equilibrium situation in the comoving frame is a static one it is characterized 
by Ui = 0. We now assume the equilibrium situation to be slightly disturbed, so that rq ^ 0. 
So we shall write 


p = po(r,z) +5p(x,y,z), p = p 0 (z) + Sp(x,y,z,t). (4) 

and treat tq, 5p , and 8p as quantities of the first order of smallness so that products of such 
quantities can be ignored. 
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We are interested in the study of the gravity wave oscillations in the disk transition layer 
where the accretion flow is sub-Keplerian and the temperature of the flow is of order 5 keV 
or more (see cartoon diagram of the system in Titarchuk, Osherovich & Kuznetsov 1999, 
hereafter TOK, Fig. 1). These temperatures of order of 5 keV are a representative values 
for plasma temperatures inferred from X-ray spectra during kHz QPO events (see more in 
Titarchuk, Bradshaw & Wood 2001 and Titarchuk & Wood 2002). 

At this stage we introduce the Cartesian coordinate system. We take the z axis to be in 
the direction of the (upward) vertical and the x axis to be such that the ( x , z) plane contains 
the angular velocity vector Q which then has components (f2 I; 0, 17*). In the general case 
the gravitational force has components g(x/R,y/R,z/R) where R = (x 2 + y 2 + z 2 ) 1 ^ 2 . 

If ( u,v,w ) are the components of the velocity u, then on substituting in equations (1), 
(2) and (3) we find 


dv 

’at 


P o 


du 

'dt 

dw 

dt 


dSp dpo < 

-w +u s +v - 


6p c 

~ C 9xVP} 

ox 

(5) 

8p 

(6) 

5p 

= -j^-g z 8p, 

(7) 

dv dw 

(8) 

dy+dz 

'o . dpo n 
-+w— = 0. 
/ dz 

(9) 


One can see that the variation of p is ignored in all terms except for the ones representing 
the buoyancy force (see C61 for details). 


3. Analysis of various types of the solutions of the problem 

In order to illustrate the buoyancy (R-T) effect we consider the simplest solutions of 
the problem including those which already exist in the literature. 

Case A 

In the case when we allow the vector fl to rotate around the 2 axis at angle <5, we 
introduce the cylindrical coordinate system, namely taking the 2 axis in the direction of the 
(upward) vertical, the x axis as the radial axis and the y axis as the azimuthal axis in the 
horizontal plane. Then the vector fi has components (Q sin 8, 0, fi cos 8) and gravitational 
force vector has components g(r/R,0,z/R ) where r = (x 2 + y 2 ) 1 ^ 2 . 
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(I) With assumptions that po is only a function of R the scale height of po is order of 
R we obtain that Sp = {(dp/dx)dx + ( dp/dz)dz ) ~ p 0 / R[{r/R)x + (z/R) C] where x = dx, 
£ = dz are the radial and vertical components of displacement respectively and conse- 
quently the radial and vertical components of the perturbation velocity are u = dx/dt , 
w = dC,/dt. Thus the vector Q = 5p(g x ,0,g z ) is transformed into Q = po{g / R)[(r / R)(rx/ R+ 
z(,/R), 0, (z/R)(rx/R + z(/R\. 

Furthermore if we neglect the effects of pressure gradients, the set of equations (5-7) 
can be reduced to the system where the right hand side consists of only the vector Q. Then, 
for radial, azimuthal and vertical displacements T, £ this system takes the form: 

X “ 2f2coe*t - —(g/R)(r/R)(r X / R + z<JR\ (10) 

T + 2f2cos(5x-2fisin<5C = 0, (11) 

£ + 2f)sin<5T = -(g/R)(z/R)(rx/R + z(/R). (12) 

OT99 have already analyzed the solution of Eqs. (10-12) (see Eqs. 2-4, in OT99) in the 
case of z/R «C 0 and 6 1. They found that in the rotational frame of reference, the 

radial oscillation with the main frequency u > k = (g/ R,) 1 '' 2 is split to the oscillations taking 
place near the horizontal (disk) plane (yT) with the hybrid frequency u>h = (w^ + 4S2 2 ) 1 / 2 
and to oscillations taking place near the vertical plane (T £) with the low branch frequency 
u>l = 2Q(u>K/oJh)sin6. 

For 2 ~ r the dispersion equation for the frequency u> 

lu 2 [u 4 — 2(u>l + 2uj 1 )lJ 1 + u* + 4f2 2 u 2 + 2fl 2 o;* sin 2 <5] = 0 (13) 

besides the nonoscillatory mode (co = 0), describes two oscillatory eigenmodes. For (5 < 1 
they are uq = u* = (g/2R) l R = u>k/V2 and w 2 = (w 2 + 4S7 2 ) 1 / 2 = {u>^/2 + 4 fl 2 ) 1 / 2 . 

The relation of the model eigenfrequencies ui = u>l/2tt, Uk = uik/2tx and v\ x — u)h/2Tr 
with the QPO observed frequencies: horizontal branch oscillation freqiencies vhbo , kHz 
frequencies were studied by OT99; TOK, Kuznetsov & Titarchuk (2002); and Titarchuk 
(2002), hereafter T02. It is worth noting that the relation between v K , Q, iq, and v L predicts 
the existence of the invariant quantity 5. (e.g. Titarchuk & Osherovich 2001). T02 calculated 
<5 and its uncertainty of 5 finding that the inferred <5— values are consistent with being constant 
at least for four Z sources, Sco X-l, GX 340+0, GX 5-1, GX 17+2 (see more on this issue 
in section §4). 

(II). With assumptions that gAp < 1 and | Vp| < 1 the set of equations for x , T, £ is similar 
to Eqs (10-12) where the right hand side vector is the zero vector, 0 = (0,0,0). In this case 
the dispersion equation for u ( for small 5 << 1) 

u> 2 (uj 2 - 4Q 2 ) = 0 


(14) 
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has the only one nontrivial root uj = 20 which is related to the eigenmode oscillations taking 
place parallel to the disk (yT)-plane. We discuss the application of this solution to the 
observation in section §4. 

Case B 

If we do not neglect the pressure gradient and assume the vertical gravitational force 
and the vector Cl directed along the vertical (8 — 0°), then the entire problem (Eqs. 5-9) is 
reduced to set of equations which have already been analyzed by Chandrasekhar in C61. He 
studied two cases: (1) In the case of two uniform fluids separated by a horizontal boundary 
he showed that for two adjacent, hydrostatic, inviscid fluids, the low fluid having density p\ 
and the upper layer having density p 2 the eigenfrequency t'R-x was 

r'R-T = [2(Q/27r) 2 + (4(fl/27r) 4 4" I/ o) 1,/2 ] 1/,2 (15) 

where is 0 is the frequency in the absence of rotation. The frequency is of order of i/k = 
(i g/R )*/ 2 if the k-wave number (i/ 0 depends on k and density difference A p = p 2 — pi) of 
order R~ l and Ap/(pi + p 2 ) ~ (0.5 - 1). 

For a stratified medium of density p = Cq exp(3z), where 1/(3 is the scale height, the 
R-T instability occurs for positive (3, whereas stable gravity waves occur for negative (3. The 
R-T frequency ^r-t = Vh = [Cj + 4(fl/27r) 2 ] 1//2 if one assumes that the wave number k is 
of order of d~ l ~ H~ l (where d is a layer size). It is worthwhile to emphasize that even 
in the simplest case of the vertical gravitational force (C61) the hybrid frequency Uh is an 
eigenfrequency (which is also true for the case with no pressure gradient effects, see above). 
Formula (10) for pr_x can also be reduced to zy, if one assumes that is 0 ~ 2(0 /2n). 

Case C 

Now we consider the case when the vector Cl rotates uniformly around the vertical at 
angle 8, then the gravitational vector Q = 8p(gr/ R,0, zg/ R), also taking into account the 
pressure effects. The general case of Q = 8p(g x ,g y ,g z ) can be analyzed in a similar way and 
it will be presented elsewhere. Following the usual practice in problems of this kind, we seek 
solutions of equations (5-9) which are of form (see e.g. C61) 

u, v, to, 8p. 8p = constant x exp (ik x x + ik y y + ik z z + nt ), (16) 

where k x , k y and k z are the horizontal and vertical wave numbers of the harmonic perturba- 
tions respectively. We also assume that po = f{x)<p(z) is a function of x (or r) and z which 
the scale heights are I/7 and \j (3 for radial and vertical density profiles respectively, r ~ R. 
Upon substituting for u,v,w,dp,8p in the form (16) equations (5) to (8) become 


nu — 2 Cl z v = —ik x (8p/ po) 4- gu'y/n + gw(3/n , 


(17) 
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(18) 

(19) 

( 20 ) 


nv — 2 Q x w + 2 Q z u = — ik y (5p/po), 
nw — 2fi x t> = —ik z (8p/p 0 ) + (zg/Rju'y/n + (zg/R)wd/n, 
k x u + k y v + k z ui = 0. 

We also use the relation nSp + (7 u + (3w)po = 0 which follows from Eqs. (9) and (16) in 
order to express 8p through po, u and w. The set of equation (17-20) assumes a nontrivial 
solution only if the determinant of the system V 0. This equation provides the dispersion 
relation for the determination of n: 

P{n ) = a 4 n 4 + a 2 n 2 + ain + ao = 0, (21) 

where a 0 = k 2 g 2 ^ 2 (k z /k x - z/R){3h ~ k z /k x ), ai = -f3k y g(2Q x k x 4- 2tt z k z )(l - 7 z/(3R), 
a 2 = 4 (Q x k x + n z k z ) 2 - g(3{-{i/(3)k 2 x {3/i - k z /k x ){k z /k x - z/R) + k 2 y { 7//? + z/R)}, a 4 = 
k 2 — kl + k 2 + k 2 z . 


3.1. Stable gravity modes 

The specific wave values k x ,k y ,k z are determined by the conditions imposed on the 
oscillatory domain boundary. Thus for a given set of boundary conditions the analysis of the 
E-T instability is reduced to the analysis of the roots of of fourth order algebraic equation 
(21), which depends on the main parameters of the atmosphere, /?, 7 and ratio of the wave 
numbers k z /k x . We assume that z/R < 1 and k x , k y oc 1/E, k z oc 1/ z is the case of interest. 

Figure 1 illustrates the specific behavior of the polynomial P(n) which should help one 
to understand the presence (or absence) of its roots for given coefficient a 0 , a 1, a 2 (o 4 > 0, 
a 3 = 0). For example, if ao < 0, ai < 0 and a 2 > 0 (see case i below) then P(n ) has two 
complex conjugate roots n\ t2 = £ ± ig and two real roots n 3 , n 4 (in fact, cPP/dn 2 > 0). 
Because ai < 0, the absolute value of the positive root n 4 is larger than that of the negative 
one n 3 . But a 3 is a sum of the polynomial (real and complex conjugate) roots, and because of 
a 3 = ni+n 2 +n3-fn 4 = 0 one can come to conclusion that 2£ = — (ri3+n 4 ) < 0. Therefore (in 
this case) the stable oscillatory mode exists and ni i2 are related to these damped oscillations. 

Below we present all cases with the oscillatory stable solutions: 

Case i: (3, (3/ 7 > 0, 3h ~ k z /k x < 0 and a 2 > 0 (see Fig. 1, curve i). For such 
conditions a 0 < 0 and a 1 < 0. Then equation (21) has two real roots which relate to the 
unstable (growing) and stable (decaying) modes with one pair of complex conjugate roots 
corresponding to the stable oscillatory mode: 


flip « a\/2d ± iuh, 


( 22 ) 
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where u>h = [(I02I + d)/2a 4 \ l C and d = (a 2 — 4a 4 ao) 1 / 2 . These roots of equation (21) are 
found using the sequential approximation method: first we solve equation (21) with ci\ = 0, 
h and then in the next stage we look for roots of equation (21) as n = h + a. 

Case ia: (3, (3 fa > 0, (3 fa — k z /k x < 0 and a 2 < 0. For such conditions ao < 0 and 
ai < 0. Equation (21) has just one pair of the complex conjugate roots and two real roots 
which relate to the unstable (growing) and stable (decaying) modes (see Fig. 1, curve ia). 
The oscillatory mode 7*1,2 is stable for which we have 

n 1,2 ai/2d± ioj L , (23) 

where wl = [( — |a 2 | + d)/2a 4 j 1 / 2 . 

Case ib: {3 > 0, (3 fa < 0, (3 fa — k z /k x < 0. For such conditions ao < 0 and a 1 < 0 
and a2 > 0. This case is similar to case i (see Fig 1, curve i) when equation (21) has two 
real roots which relate to the unstable (growing) and stable (decaying) modes and complex 
conjugate roots which correspond to the stable oscillatory mode: 

7*1,2 ~ a\/2d i ioJh- (24) 

Case ii: f3 > 0, (3fa — k z /k x > 0. For such conditions a 0 > 0 and < 0 and a 2 > 0. 
Equation (21) has a pair of complex conjugate roots (see Fig. 1, curve ii) and one of them 
7*1,2 is related to the stable oscillatory mode: 

711,2 ~ ai/2d ± iuJh- (25) 

Case in: (3 < 0, (3fa — k z jk x > 0, and 02, d > 0. For such conditions ao > 0 and aj > 0. 
Equation (21) has a pair of complex conjugate roots (see Fig. 1, curve iii) and one of them 
713,4 is related to the stable oscillatory mode: 

m, 2 « —a\j2d ± iujL, (26) 

where lo l = [(a 2 - o?)/2a 4 ] 1 / 2 . 

Finally, we single out the case when the effective gravitational force goes to zero. This 
can happen during a burst event, when the gravitational forces are compensated by the 
radiation pressure forces (e.g. Titarchuk 1994). For such conditions, ao, ai — ► 0 and a 2 /a 4 = 
4[(A : x /k)Q x /k+ (k z /k)Q z ] 2 . Equation (21) has zero roots n = 0 and conjugate complex roots 

7ii.2 = ±i(a 2 /a 4 ) l/2 = ±2 i[(k x /k)Q x + (k z /k)£l z ] (27) 

which are related to a pure harmonic mode. In fact, this result also follows from treatments 
detailed in C61 and OT99. The hybrid frequency then becomes simply v h = 2(fl/27r) for 
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g = 0 . This case of the perfect coherent oscillations would also occur if the density profile 
is quasi-uniform i.e. when /?, 7 — > 0 but p/'y = 0 ( 1 ). These two conditions (on either the 
effective gravity force or the density profile) can be realized during the burst event. 

In this section the main goal is to reveal all cases when stable gravity modes exist 
and when the stability breaks down. This analysis is particularly important in a view of the 
transient nature of QPO features (see e.g. Zhang et al. 1998 ). With an increase of bolometric 
luminosity (presumably in mass accretion rate) the kHz QPO frequency increases and then 
entirely disappears! At low rates it appears once again! The stability analysis presented 
here leads to conditions for the existence and the destruction of gravity modes in terms of 
density profiles scale heights 0 ~ l , 7 _1 and boundary conditions (k x , k y k z ). In fact, the 
strong dependence of the gravity wave stability on the density profile was a central point of 
Chandrasekhar’s analysis (C 61 ). For the accretion disk cases, we can give an example where 
stable gravity modes are followed by instabilities (the gravity mode destruction). In case 
(ib) we have a stable mode with tOh as a QPO frequency. When the density profile changes, 
over 2 — coordinate from 7 < 0 to 7 > 0 then the stable g-mode with uJh can still be sustaned 
[see case (ia)], But if the density profile over z coordinate is stabilized (ft < 0 ) and that over 
r coordinate is inverted (i.e. from 7 < 0 to 7 > 0) then the QPO oscillations are no longer 
stable. 


4. Gravity modes and their relation to QPO phenomenon 

As we have seen in the previous section, there are two stable oscillatory gravity modes: 
one is associated with the hybrid frequency ujh (see cases i, ib, ii ) and another with the low 
branch frequency ul (cases ia and in). We also estimate the decay rate for oscillations as 
A = \a x \/d and the QPO quality value Q = u>j 2A. The presence and absence of these modes 
depend on the atmospheric structure (scale height inverses (3 and 7) and on the imposed 
boundary conditions (wave numbers k x , k y and k z ). Furthermore, it is possible to restore the 
related boundary conditions if one compares the observed QPO features with the calculated 
mode frequencies and Q- values. The analysis made in section 3 is also necessary in order to 
control an accuracy of numerical calculations of the set of hydrodynamical equations ( 5 - 9 ). 

To illustrate the results obtained in §3 we should specify the orders of the introduced 
quantities k x , k y , k z ,/ 3 , 7. Namely we assume that k x ~ 1 / R, k y ~ 1 / R, k z ~ 1/2,7 ~ q/R 
(where |g| < 1) and 2 / R < 1. If we also assume that 0 /y —k z /k x ~ R/z then 0 ~ 2< 7/2. With 
these assumptions expressions for the polynomial coefficients aj (see Eq. 21 ) are significantly 
simplified: a 0 /a 4 ~ q 2 uj^, a x /a 4 ~ — 2 \q\(Jy i {zVt x / R + fl z ), a^ja^ ~ 4 (zQ, x /R + fi 2 ) 2 + qu^. 
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Because |(/| < 1 and lok < Uh we can calculate as 

~ I^K^k/^/iVk (28) 

and 

Q l = uj L /2(a l /2d) « u h /2Q > 1/2. (29) 

TOK introduced the classification scheme for the QPO features and related the observed 
high and low kHz frequencies to the hybrid and Keplerian frequencies, = u>h/ 2ir and 
r'K = u>k/2tt respectively. They also attributed the hectohertz frequencies detected in the 
atoll source 4U 1728-34 (Ford & van der Klis 1999) to the low branch. The angle 5 was found 
to be almost twice that found in Sco X-l (see also T02 for details of 5— determination). The 
low frequencies in atoll sources (as a rule) are three times higher than those in the Z-sources. 
The hectohertz frequencies have been identified in several other neutron star LMXBs (4U 
0614+09; van Straaten et al. 2000, 2002; SAX J1808.4-3658, 4U 1705-44; Wijnands &; van der 
Klis 1998). They weakly depend on the kHz frequencies, their ratio to low kHz frequency 

being about (4-5) (van Straaten et al. 2002, hereafter S02). The hectohertz frequency 

Lorentzian profile are very broad with (/-values around 1 or even less. We suggest that these 
observed frequencies can be identified as the low branch frequencies (see formulas 28-29). 
Taking into account resonance effects (T02), we correct the error bars of the hectohertz 
frequencies presented in S02. In Figure 2 the best fit to the data (S02) is presented using 
formula (28) which includes one fit parameter q. The best-fit value of jqr| = 0.3 for which 
Xred = 0-76. The estimated QZ 1/2 are also in agreement with the observed values of Q (S02, 
Table 2). 

The harmonic modes u>h m = 2(k x Q x + k z Q z ) can be related to a coherent oscillation 
frequency of 582 Hz which is observed in 4U 1636-53 during the superburst (Strohmayer & 
Markwardt 2002). If k z /k x 1^ 1 (which is our case) is a double frequency of the NS spin 
frequency. Furthermore, we have already found the same eigenmode with the frequency 2S1 in 
case (All) when effects of the buoynant forces and pressure gradients are neglected (see §3). 
In fact, Miller (1999) has reported the detection of the coherent pulsations of frequency ~ 291 
Hz during the burst development in 4U 1636-53. Thus the NS spin frequency Q/2ir = 291 
Hz and the double NS spin frequency 2Q,/2 tt = 582 Hz have probably been detected in 4U 
1636-53. 

It is also worth noting that the RXTE observations of neutron binaries establish that 
the kHz QPO frequencies do not seem to exceed beyond a certain upper limit (Zhang et al. 
1998). This observational effect may be a result of the density profile. Our stability analysis 
presented in §3 clearly indicates such a possibility (see also C61). 
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5. Conclusions 

I have presented a detailed study of the Rayleigh-Taylor (R-T) instability in the accre- 
tion flow. To summarize I have :(1) put forth arguments to explain the QPO phenomena, as 
a result of the R-T effect in the rotational frame of reference. (2) formulated and solved the 
mathematical problem of the gravity wave propagation (R-T effect) in the accretion flow. 
(3) concluded that the stable gravity modes in the rotational frame of reference are related 
to the hybrid and low branch frequencies. (4) demonstrated that the particular problem of 
the gravity waves in the rotational frame of reference, in the approximation of very small 
pressure gradients, is reduced to the problem of the classical oscillator in the rotational frame 
of reference which was previously introduced and applied for the interpretation of kHz QPO 
observation by OT99. 

I demonstrate that these frequencies are intrinsic features of the R-T effect. They ap- 
pear in various configurations of the accretion flow depending on assumptions regarding the 
density profiles, the boundary conditions and the effects of the pressure forces. It is not by 
chance the high and low frequencies phenomenon has common observational appearances for 
a wide range of objects classes from black hole sources down to white dwarfs (Mauche 2002). 
(5) Investigated the conditions for the density profile and the wave numbers (boundary con- 
ditions) when the gravity modes are stable. (6) Identified the observed QPO frequencies 
seen in the power density spectra of NS LMXBs using the inferred gravity mode frequencies. 
In particular, I found that the inferred low branch frequencies and their Q-values are consis- 
tent with the QPO hectohertz frequencies observed in the atoll sources 4U 1728-34 and 4U 
0614+19. (7) During the NS long (super) burst event, I find that the observer should see 
oscillations at double NS spin frequency. The Coriolis force is the only force which acts in 
the rotational frame of references and its presence causes perfect coherent pulsations with a 
frequency twice of the NS spin frequency. 

Finally one can conclude that the R-T gravity wave oscullations must be present and the 
related QPOs should be detected in any system where the gravity, buoynancy and Coriolis 
force effects cannot be excluded (even in the Earth and solar environments). 

L.T. acknowledges fruitful discussions with Chris Shrader and Kent Wood. I also ap- 
preciate the fruitful discussions with the referee and his/her constructive evaluation of the 
manuscript. 
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Fig. 1. — Behavior of the dispersion polynomial P(n) for stable gravity modes (see text). 
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Fig. 2. — Correlation between the lower kHz frequency and hectohertz frequency for 4U 
1728-34 (star), 4U 0614+09 (triangles) (van Straaten et al. 2002) and the best-fit curve of 
low branch frequency vs Keplerian frequency (solid line). Xred. = 0-?6 for \q\ = 0.3. 




